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This prehminary report proposes integrating the Maxwell equations in Minkowski spacetime us- 
ing coordinates where the spacelike surfaces are hyperboloids asymptotic to null cones at spatial 
infinity. The space coordinates are chosen so that Scri-|- occurs at a finite coordinate and a smooth 
extension beyond Scri+ is obtained. The question addressed is whether a Cauchy evolution numer- 
ical integration program can be easily modified to compute this evolution. In the spirit of the von 
Neumann and Richtmyer artificial viscosity which thickens a shock by many orders of magnitude 
to facilitate numerical simulation, I propose artificial cosmology to thicken null infinity Scri-|- to 
approximate it by a de Sitter cosmological horizon where, in conformally compactified presentation, 
it provides a shell of purely outgoing null cones where asymptotic waves can be read off as data on 
a spacelike pure outflow outer boundary. This should be simpler than finding Scri-|- as an isolated 
nuh boundary or imposing outgoing wave conditions at a timelike boundary at finite radius. 



It is a pleasure to dedicate this work to my friend and 
colleague Stanley Deser as we celebrate his ancient works 
and his continuing contributions to a broad variety of 
problems in physics. 



I. INTRODUCTION 

As LIGO 1] approaches its first stage design sensitivity 
and other |2i gravitational wave observatories progress, 
theoretical descriptions of the expected waveforms from 
the inspiral of binary black holes and neutron stars re- 
main less detailed than one might hope. The largest ef- 
forts are devoted to solving Einstein's equations by dis- 
crete numerical methods based on a space-|-time split of 
the equations and intend to impose as a boundary con- 
dition an asymptotic Minkowski metric. Initial formu- 
lations were given the ADM name although (probably 
inconsequentially) none uses the conjugate sets of fields 
gij and tt*-' favored by ADM. Variations on these formu- 
lations have shown practical improvements |^, ^ HI , and 
others Q with less developed implementations hold the 
theoretical suggestion of better numerical behavior. 

Because the (expected) observable gravitational waves 
will be seen far from their sources, some efforts have asked 
that the data set representing the state of the gravita- 
tional field at any computational step be not that on an 
asypmtotically flat spacelike hypersurface, but instead 
that on an outgoing light cone. Thus the Pittsburgh 
group j2i| has focussed on slicing spacetime along null 
cones. By extracting a conformal factor as in Penrose 
diagrams, null infinity becomes a finite boundary in the 
computational grid. But these outgoing null cones do 
not have a useful behavior near the wave sources, lead- 
ing to a study matching them to flatter spacelike slices 
there. Friedrich jSj] formulated the Einstein equations in 
a way where the conformal factor could be one of the dy- 



namical variables, and also considered the case where the 
time slices were spacelike but asymptotically null, as for 
hyperboloids in Minkowski spacetime. A large program 
to implement this approach is being pursued by Sascha 
Husa jO, 10] and others at AEI-Potsdam. 

This paper aims to formulate within the simpler prob- 
lem of solving Maxwell's equations some aspects of these 
methods which could be attacked with smaller computa- 
tional resources. In particular we consider the question 
of whether fairly straightforward evolution algorithms de- 
signed for spacelike slicings of spacetime will encounter 
special problems when those slices are asymptotically 
null, and the question of imposing boundary conditions 
at null infinity (Scri+). 

Work with Mark Scheel and Lee Lindblom con- 
sidered the first of these questions, but employed the 
scalar wave equation. Those preliminary results were 
presented at the KITP Workshop "Gravity03" in June 
2003. Vince Moncrief ^3 then suggested that since the 
coordinate choice which brought Scri-|- to a finite radius 
displayed a conformal metric which was nonsingular at 
Scri-|-, it would be better to study instead the Maxwell 
equations whose properties under conformal transforma- 
tions are simpler that those of the scalar wave equation. 

A new suggestion is given below that the boundary 
conditions at Scri-f might be easier to handle if Scri+ 
were approximated by a de Sitter cosmological horizon 
with a cosmological constant chosen for numerical con- 
venience rather than for physical interest. 



II. MAXWELL'S EQUATIONS 

In formulating a 3-1-1 statement of the Maxwell equa- 
tions I, of course, revert to lessons from Arnowitt and 
Deser [I3. Following their Schwinger tradition I use a 
first order variational principle SI = with 
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where 



and 



ui/a/3 _ 1 

y — 2 



(2) 



are components of an operator that maps covariant anti- 
symmetric tensors onto contravariant anti-symmetric 
tensor densities. This mapping has an inverse whose com- 
ponents are 



gt^iyap = \{gti.agvp - gpipgva)/ V-g 



(3) 



Note that these two symbols are not related merely by the 
usual lowering of indices, as one map raises the density 
while the other lowers it. Their product 



(4) 



gives the identity mapping of antisymmetric tensors onto 
themselves. 

In this variational principle the fields and F^^i, are 
varied independently and yield the equations 



and 



Here 



-gF'^" arises as 



(5) 
(6) 
(7) 



Again following the given tradition I seek a proto- 
Hamiltonian form pq — H for this Lagrangian and rewrite 
it by replacing the field in / by its equivalent 5^^" 
field, giving 



(8) 



The only term here containing time derivatives is 
-S^'doAi so that 



(9) 



and Ai are the only fields for which one obtains evolution 
equations. The scalar potential 



i' = ~AQ 



(10) 



and the are therefore Lagrange multipliers which en- 
force the constraints 



as* = 



and 



(11a) 
(lib) 



-^ij — gij,af!-S°'^ — diAj — djAi 

The evolution equations can then be written as 

doA, = -E, - (12a) 



where 



daD' = djSj'' (12b) 



= 5,o,a;3r'' and ^'^=r' ■ (13) 



There is no evolution equation for ijj = ~Aq but one 
can be supplied as a gauge condition. There is also no 
evolution equation directly from the variational principle 
for Bij but one can be deduced by taking the time deriva- 
tive of equation l|llb|l and evaluating the time derivatives 
on the right hand side using equation (|12a|l . The result 



IS 



'do Bij — diEj — djEi 



(14) 



In a numerical example mentioned below, the fields 
were taken to be Ai and using the evolution equa- 
tions H12a() and (jl2b|l . Equation ljllb|) was treated not 
as a constraint but as merely making Bij an abbreviation 
for the right member of that equation, thus introducing 
second (spatial) derivatives into the system of equations. 
For that example of a wave packet the gauge chosen was 
Coulomb gauge with ^ = for all time. The remaining 
constraint (|lla|l is easily seen to be preserved in time as 
a consequence of the evolution equation (|12b|l . 

Another possible system of equations would be to use 
(jl2b|l and (|14|l as evolution equations. Since the vector 
potential does not appear in these equations, it can be 
ignored and not evolved. The constraint (jllbp would 
then have to be replaced by its integrability condition 



d[kBij] — 



(15) 



treated as a constraint equation. It is easily seen from 
(|14|l that this constraint is preserved in time (differen- 
tiably) as a consequence of (|14|l . The Gauss constraint 
(|lla|l remains part of this system and is preserved in time 
by exact solutions. These equations l|12bl) and 1)14(1 are 
exactly what one writes in classical electromagnetism in a 
material medium (e.g., 0, Section 16-2]); only the con- 
stitutive equations (|17|l are new. With the usual 34-1 
decomposition of the metric 



these are 



-a^dt^ + jijidx' + I3'du){dx^ + (i^ du) (16) 



and 



(17a) 



(17b) 



which are more complex than the E = D/e and H = 
B//Lt of simple isotropic media. The Einstein ^ther is, in 
coordinate terms, anisotropic not only in space but also 
in spacetime, but conveniently linear. 
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III. CONFORMAL INVARIANCE 



IV. BACKGROUND SPACETIME 



This exploration of numerical evolution with the 
Maxwell equations was provoked by Vince Moncrief's 
reminder that the Maxwell and Yang-Mills equations 
(in four spacetime dimensions) have the simplest pos- 
sible conformal structure. A confornial transformation 
g^i, 1-^ ^"^g^v changes the metric geometry to a differ- 
ent one which shares the same light cones. Under this 
change the metric dependent factor g^^^'^l^ in the vari- 
ational integral / does not change since g^'' ^ fl^^g^"^ 
while ^/—g i— *■ fl'^y^—g. As nothing in the variational inte- 
gral changes under a conformal transformation, it follows 
that if the fields and i^^^, or ^^"^ satisfy the field equa- 
tions in one metric, so do they also in any conformally 
related metric. 

Although nothing in section ^1 depends on the metric 
having any special properties (other than Lorentz signa- 
ture) such as flatness or satisfying the Einstein equations, 
these conformal properties in four spacetime dimension 
are particularly useful since the metric H21I) we intend to 
use is singular only in a conformal factor s^/g^ which we 
now see will appear nowhere in our field equations. A 
similar simplification occurs when the metric of interest 
is the Schwarzschild metric in the hyperboloidal slicings 
used in 

In the 3-1-1 decomposition of the field equations we note 
that AijT)'^, Eijtp, Bij , i^'^ and are conformally invari- 
ant, while 7y I— > f2^7y, a i— > fla, i— > fl^^/j. Thus the 
metric structures which appear in the constituitive equa- 
tions ((T7|l . {a/^)jij and 0^7**^7^' as well as are 
conformally invariant. 

The constraint equations Hll|l and (|15|l have an even 
stronger invariance: they are metric invariant. No met- 
ric quantities appear in these constraint equations when 
they are written in terms of our choice of fields. Thus 
if two hypersurfaces are diffeomorphic then initial con- 
ditions for the Maxwell equations which satisfy the con- 
straints in one Lorentzian manifold can be imported to 
the other where they will again satisfy the constraints. If 
the two manifolds are not conformally equivalent, how- 
eve, the subsequent evolutions of these initial conditions 
will generally be inequivalent. 

Gauge conditions present greater difficulties. In the 
absence of sources the "Coulomb" gauge ip = is con- 
formally invariant, but only 3-dimensionally coordinate 
invariant. On the other hand the simplest Lorentz gauge 
d/j.iy^—gg'^'' A^) = is coordinate invariant but not con- 
formally invariant. The set of field equations which 
evolve J)* and Bij without the use of the vector poten- 
tial are conformally invariant but do not require a gauge 
condition. 



The background metric for this study will be the de 
Sitter spacetime 

ds^ = -dT^ + dX^ + ^ ^^2 _^ (^ji^ I L^)[cLT -dRf . 

(18) 

which is here presented in Kerr-Schild form as fiat space- 
time plus the square of a null form. The cosmological 
constant here is A = +2>/L'^. Since dT — dR is null in 
both the underlying flat spacetime and after the Kerr- 
Schild modification, the hyperboloidal slices we used in 
fiat spacetime should remain asymptotically null also in 
this modified metric. Therefore we introduce the same 
coordinate change as in our earlier scalar work 



and 



T 
s 

s 



1„2 
2' 



It-2 
4' 



i„2 
4' 



where = X^+Y^+Z^ = X'X' and = 
x'x*. This leads to the metric 



ds' 



where 



isyq^)[-a^du^ 



(19) 



(20) 



(21a) 



L2(i 



{dr + Pduf + r^d^^] 



g = 1 



1^2 



L2(l + ir)4 



and 



L2(H-ir)2 



1 



i2(l-f ir)4 



(21b) 



(21c) 



(21d) 



I have verified that this metric (|21|l actually is the de 
Sitter metric by two tests: I have run a GRtensor II 
Maple 8 worksheet to calculate the Einstein equations for 
it and find 



(22) 



and I have also had that worksheet calculate the Riemann 
tensor which is 



R^"^p = {i/L^Wj'{,-5';5';:) 



(23) 



and shows that this is a spacetime of constant curvature, 
therefore de Sitter. Richard Woodard [3| has confirmed 
this by showing that the metric H18|) can be reduced to the 
familiar static deSitter metric by the coordinate change 
T = f + f{R) where df/dR = -R^/{L^ - i?^). 
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V. CAUSAL STRUCTURE 

This de Sitter metric is presented in equations 1)21(1 
with an apparent singularity at r = 2 where g = 0. This 
is probably a coordinate singularity since the curvature 
does not go bad there, and also because the de Sitter 
spacetime in notorious for being presented in many differ- 
ent guises corresponding to different coordinate systems, 
many of which cover only a patch of the full manifold. 
But for present purposes we only need this coordinate 
patch, and the singularity at i? = oo will be removed by 
a conformal transformation later. 

By dropping the conformal factor in the metric 121() 
we arrive at a nonsingular metric 



-a du' 
- 1 + 



(24) 



i2(l + ir)4 



{dr + Pduf + r^d^^ 



using the same abbreviations as in equations H21|l . Note, 
however, that a and (3 are now the lapse and shift of 
this regulated metric. This metric no longer satisfies the 
Einstein equations (with A term) but is related to the 
original (de Sitter) metric in a known way. From this 
point on we use only this metric, and the previous met- 
ric will be called ds^ = fi^ds^ following the practice in 
Wald's discussion I?', Appendix D] of conformal trans- 
formations, with ^ — s/q in our case. 

The causal structure of the regulated metric (|24|l is 
simpler than that of the original metric (|18|l . The hyper- 
surfaces of constant u are everywhere spacelike since the 
regulated metric is positive definite when du = 0. The 
hypersurfaces of constant r are spacelike only in an in- 
terval around r = 2 as seen from their normal vector dr 
being timelike there: 



(Vr) 



(1 



1,2)2 



{rs/LY 



(1 



l^2^2 



(25) 



This norm is plotted in Figure^ over a broad range and 
again in Figure El giving the interval around r = 2 in 
more detail. 



VI. ANALYSIS 

What can we learn using the tools described above? A 
few lessons already appeared in the work 11] with Scheel 
and Lindblom. Suppose only the time coordinate were 
changed to parameterize a family of hyperboloids, here 



[T - S{U - 1)]2 - i?2 = s2 



(26) 



in the original coordinates of equation H18|l . Then the 
outgoing coordinate speed of fight dR/du would become 
unboundedly large as i? — > oo. The Courant condition 
usually met in numerical implementations would then 
make the allowable time step impractically small. This 
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FIG. 1: The squared norm (Vr)^ = g^'~ of the normal Vr 
to hypersurfaces of constant r in the conformally regulated 
metric 1241 is plotted over a range from the origin to beyond 
the location r = 2 of the surface which becomes Scri-|- in the 
high resolution limit L — > oo. This example uses L/s = 10. 




FIG. 2: The squared norm (Vr)^ — g^'~ of the normal Vr to 
hypersurfaces of constant r in the conformally regulated met- 
ric H24|l is plotted over a narrow range arround r = 2 where 
r = const are spacelike pure-outflow hypersurfaces on which 
the entire forward light cone is oriented toward increasing r 
when g^'^ < 0. This example uses L/s = 10. 



can be cured by "Analytic Mesh Refinement" (AnMR), 
which is just a coordinate transformation such as that 
from i? to r in equation (|20|1 . since outgoing waves get 
redshifted as the u = const hypersurfaces tend toward 
outgoing null cones. Because this redshift gives long 
wavelengths, high resolution is not needed and the large 
R domain can be compressed to a finite range. Thus we 
find that the outgoing coordinate speed of light dr/du 
becomes 



Cout 



(27) 



which remains less than 4 out to r = 2 or i? = oo. The 
ingoing radial light speed Cin is similarly bounded in mag- 
nitude but becomes positive in the region where (Vr)^ is 
negative rather than the simpler value Cin = —(I — ^t)^ 
in the Minkowski (L = oo) case where it merely touches 
zero on Scri-|-. Thus the encouraging results from 
are unchanged by introducing the cosmological term. 
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VII. EXACT SOLUTIONS 

There is in principle no difficulty in finding exact so- 
lutions to these equations in Minkowski spacetime. One 
merely takes solutions from any textbook and performs 
the coordinate transformation to our uxyz compactified 
hyperboloidal slicing coordinates. In detail is it a great 
convenience to have GRTensor 15] available to perform 
some of the algebra and much essential checking. 

Exact solutions play two roles. The first is that they 
provide initial data for numerical integrations whose 
physical interpretation is known. Secondly, though com- 
parison, they allow testing the accuracy of the numerical 
solution. Indeed, the whole point of carrying out these 
numerical evolutions is to explore in a simpler context 
than full general relativity different approaches for pos- 
sible adaption to Einstein's equations. Thus one may 
explore and compare different integration schemes either 
to understand them better as in the work of Baumgarte 
[Tsj l , or to test possible new schemes as in . Questions 
which may be treated in this simplified context include 
the stability of different formulations, the control of con- 
straints, the treatment of boundary conditions, the ex- 
traction of asymptotic wave amplitudes and waveforms, 
etc. Here we address the use of hyperboloidal time slices 
and the possibilities for carrying the integration out to 
and beyond Scri-|- using a finite grid. 

We find that a sufficient variety of examples can be 
found by asuming that the vector potential A — A^^dx^ 
has the form 

A = f sin^ ed<t)~^dT . (28) 

Three cases are (1) a uniform magnetic field with ■;/; — 
0, / = BzR^ (2) a static dipole electric field with / = 
0, V' = Z/R^ (3) the Baumgarte [H choice of a wave 
packet, which is our principle test case, with ijj — Q and 

/ = (^-2AC/) exp(-A;72)-(l+2AF) exp(-AF2) (29) 
It K 

where U and V were defined by 

U = T-R = s(^u- Yq~T~) (^0^) 



and 

V = T + R = s(u+ Y~i~) ■ '^^^'^^ 

In spite of first impressions one finds that / = 0(i?^) at 
small R or small r, but numerical implementations must 
avoid evaluating this / at r = as subtle cancellations 
occur. To convert equation l|28|l to uxyz coordinates re- 
quires 

-dT ^du+ j-^r^dx' (31) 

S (1 _ 1^2)2 

and, from X' — sx'^/{l — \r'^) and the standard rectan- 
gular to spherical coordinate transformation, 

sin^ed(t) ^ X^iXdY ~Y dX) = ^(xdy-ydx) . (32) 

The results are nontrivial even for a constant magnetic 
field in these coordinates which can be useful for debug- 
ging. To obtain initial conditions for numerical integra- 
tions requires that from these vector potentials one com- 
pute the fields F^,^, and thus = Figure El 
plots an analytic solution obtained in this way at three 
different times. The incoming wave packet (at u = —2) 
is blueshifted relative to the normals the to constant u 
hyperboloid, while the outgoing packet (at u — 1.4) is 
redshifted so that high spatial resolution is not required 
to describe it near Scri-|- (which is at r = 2). Note that 
there are no irregularities of the outgoing packet as it 
crosses Scri-|-. 
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FIG. 3: For the Baumgarte choice of EM wave packet from equations 1281 1 and 12911 we plot the conformally invariant electric 
field along the i-axis at three different hyperboloidal times u — —2,0, and 1.4 . Note that the field values are finite and 
generally nonzero at Scri+ which is x — 2. Those at the latest time have been here magnified by a factor of 10. 



O.St 




-0.2 J 
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as in Figure 13 we plot ©"as function of it at Scri+ {x — 
2). Note that at constant r one has T = sit + const so this 
waveform is not distorted from its (T, R) presentation as is 
the pulse shape (Figure when plotted as function of r. 
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